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Abstract 

Fluid mechanics can be formulated on dynamical surfaces of arbitrary co-dimension embedded in a 
background space-time. This has been the main object of study of the blackfold approach in which 
the emphasis has primarily been on stationary fluid configurations. Motivated by this approach we 
show under certain conditions that a given stationary fluid configuration living on a dynamical surface of 
vanishing thickness and satisfying locally the first law of thermodynamics will behave like an elastic brane 
when the surface is subject to small deformations. These results, which are independent of the number 
of space-time dimensions and of the fluid arising from a gravitational dual, reveal the (electro)elastic 
character of (charged) black branes when considering extrinsic perturbations. 
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1 Introduction 

Recently, the interest in the study of fluid dynamics has increased considerably from the point of view 
of gravitational physics and holography. This has been triggered by different approaches to phenomena 
encountered in black hole physics leading to several long- wavelength effective descriptions. The membrane 
paradigm describes the behavior of black hole horizons when subject to certain perturbations [H El El HI E] . 
In the context of AdS/CFT, the fluid/gravity correspondence O [7| provides an effective description for 
the behavior of perturbations of black branes along boundary directions. The blackfold approach [SlI^ITUl 
[TH I12j. an effective theory for the dynamics of thin black branes, characterizes the behavior of strained 
black brane geometries along worldvolume directions and orthogonal to these All these approaches show 
that the physical properties of black holes in certain regimes can be captured by the physics of fluid flows 
that either live on the stretched horizon (membrane paradigm), on the boundary (fluid/gravity) or in 
an intermediate region (blackfold approach). These insights motivate the search for a generalization of 
the mechanics of fluid flows that live on an arbitrary surface of any co-dimension placed in a background 
space-time. 

This development has taken place within the realm of the blackfold approach |10| [T3l [H] , and has 
lead to a vast study of large classes of fluids with particular focus on stationary configurations. To make 
the reader familiar with this effective theory we present the equations of motion that describe it, first 
derived by Carter [T3] . 



Here Da is the covariant derivate with respect to the induced metric 'jab of the surface on which the 
fluid lives while Kab^ is its extrinsic curvature. When the stress-energy tensor T"'' is assumed to be 
of the perfect fluid form, the first set of equations give rise to the usual energy density continuity and 
Euler equations of a perfect fluid. When considering, for example, the type of (intrinsic) hydrodynamic 
perturbations of black hole horizons encountered in the context of the fluid/gravity correspondence, this 
set of equations can be derived as constraint equations directly from Einstein equations [6]. One can 



T^'Kab' = . 



(1.1) 
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think of perturbations of this type as fluctuations in the fields that characterize the material - in this case 
the fluid - that lives on a space-time surface. Proceeding order-by-order in perturbation theory results 



in dissipative corrections to the stress-energy tensor. The second set of equations in (1.1) is of extrinsic 
nature and is associated with deformations of the geometry (surface) on which the fluid flows. As in the 
case of hydrodynamic fluctuations, this set is also directly derivable as constraint equations from Einstein 
equations [Sj, il2j when deformations, for example bending, of black brane geometries are considered. Our 
aim in this note is to show that this extrinsic set of equations can be viewed as a relativistic generalization 
of elasticity theory of thin branes. 

The reader may be familiar with Eqs. (1.1) when dealing with Dirac branes where T*^^ = T^)^^""^ with 
being the tension of the p-brane but we note that they have a direct non-relativistic analog when 
considering deformations of thin membranes. To be precise, suppose that we are given a thin elastic 
membrane of thickness ro and subject it to external forces applied at its circumference causing it to 
stretch in all directions. Assuming the material to behave elastically, generating internal stresses a""^ 
that encode its Hookean response to the stretching, the equations of motion that govern its mechanical 
equilibrium can be obtained by varying the free energy [16] 



FfX^^] = \ I dVa'^^Uab , (1.2) 
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Vol 



where X'^ is the set of mapping functions that parametrize the position of the membrane in the ambient 



space while U""^ is the strain tensor. The resulting set of equations is exactly (1.1) with T""^ replaced 
by a""^ and with the indices a, b only running through the spatial directions. Within this perspective, a 
response of the Dirac brane type can be seen as the isotropic stretching of a p-brane. 

The motivation for this interpretation has its origin in the analysis of the corrections associated with 
extrinsic perturbations of black branes performed in \n\ \T2\ [T8] which takes into account the finite 
thickness of the brane geometry. In fact, for the particular case of black strings bent into circular 



shape, Eqs. (1.1) keep the same form, with the stress-energy tensor acquiring a correction that can be 
interpreted as a Hookean elastic response of the string due to the bending |17l I18j . In line with this, 
we ought to show that the extrinsic dynamics is already to leading order that of an elastic solid. Even 
though greatly inspired by black brane physics the scope of the results that will be presented here is far 
wider. As a matter of a fact, what will be shown is this: given a stationary fluid configuration satisfying 
local thermodynamic laws living on a thin space-time surface, assuming diffeomorphism invariance and 
ignoring backreaction, its extrinsic dynamics in transverse directions to the surface correspond to that of 
an elastic brane. We expect this interpretation to lead in the future to a formal development of general 
relativistic elasticity theory of branes that can be captured from gravity in the same spirit as a rigorous 
development of fluid and superfluid dynamics was captured using the fluid/gravity correspondence (see 
[m 1201 mi [H [221 [23]). Refs. dZl [m [ig already made important steps towards this goal by measuring 
for the first time response coefficients, such as elastic and piezoelectric moduli, of materials - black branes 
- bent both in space as well as in time. 

This note is organized as follows. In Sec. [2] we review and reinterpret the equations of motion for sta- 
tionary fluids living on dynamical surfaces [9, 10l[T3l[T3] in terms of elastic solid concepts. We show this 
to be the case for isotropic fluids that in addition may carry different types of charges. This solid interpre- 
tation along extrinsic (transverse) directions to the embedding surface also holds for the anisotropic fluids 
of \13\ I14j carrying string charge. In Sec. [s] we give examples of fluid configurations arising from gravita- 
tional duals satisfying these properties and apply the formalism with this new interpretation obtaining 
the elastic equilibrium condition for certain black holes. In Sec. |4]we conclude. 
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2 Stationary fluids on dynamical surfaces 



In this section we review and reinterpret the theory describing stationary fluid configurations confined 
to a (p + l)-dimensional dynamical surface, parametrized by the mapping functions embedded in an 
ambient D-dimensional space-time with metric g^u O E] ■ When the surface and the space-time are of 
equal dimension we obtain the ordinary description of fluid mechanics. We work under the assumption 
that the fluid is in local thermodynamic equilibrium. This is achieved when the mean free path charac- 
terized by the inverse of the local temperature 7~(cr'*) is much smaller than the radius of curvature of the 
embedding geometry i?(cj"), i.e., 

1 



where a"", a = 0,...,p are the coordinates that parametrize the worldvolume Wp+i traced out by the 
surface in space-time and endowed with metric 'jab = OfiudaX'^dhX'^ , i/ = 0, D — 1, with Lorentzian 
signature. We further assume that the surface is infinitely thin. Finite thickness effects have been 



considered in [T71 [18]. However we note that as in the classical example discussed around Eq. (1.2) 
where the thickness can be integrated out and the object replaced by an infinitely thin elastic membrane, 
the thickness of the {p + l)-dimensional objects considered here can also be integrated out giving rise 



to a worldvolume effective theory. Sees. 2.1 and 2.2 review this theory and the intrinsic fluid dynamics 



following Refs. [9lll0 t [T ^ ll3j. while the remaining reinterpret the extrinsic dynamics in terms of relativistic 
elasticity theory. 

2.1 The worldvolume effective theory 

Assuming the fluid not to backreact onto the background, the usual equations of fluid mechanics can be 
derived by imposing conservation of the stress-energy tensor T^'^ [15j 

y^ffiu ^ jc-M ^ ^2.2) 
where we have included the possibility of an external force J^^. The dynamics of fluids living on dynamical 



surfaces also follows from (2.2) but the stress-energy tensor characterizing these configurations is confined 



to the surface in the following way: 



T'"'{x'^) = / dP+V7=7r^^(a") ' ^ " . (2.3) 

Here one should interpret T^'^{x°') as the full spacetime stress-tensor while the components T^'^{a"') with 
support on Wp+i should be seen as worldvolume densities of stress-energy. We assume the force term 



J^^{x°') to be of the same form as (2.3) with appropriate worldvolume density J-^{a"'). Introducing 



2.3[) into Eq. (2.2) results in a worldvolume theory where the densities T^'^{(j°') only have components 



tangential to the worldvolume 

1_P^T^''' = , (2.4) 
where the orthogonal projector to the worldvolume is defined as 1_p^= gp'^ — jP'^ with 7^^^ = j'^'^daX^dhX'^ 



being the push forward of the induced metric. The equations of motion derivable from Eq. (2.2) can be 
written in the form 

Yu^pT^" = , (2.5) 
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together with the boundary condition 



p+1 







(2.6) 



where fii, 

on T'^'^ requires T^' 



is a unit normal vector orthogonal to the worldvolume boundary. Imposing the constraint (2.4) 
T°'^UauY with Ua = daX^ which can then be used to split Eq. (2.5) into two sets 



of equations by projecting along tangential and orthogonal directions to the worldvolume: 



ab 



(2.7) 



These equations reduce to Eqs. (1.1) when = 0. The first equation expresses the conservation of 



worldvolume stress-energy while the second can be interpreted as the balance of forces acting on the fluid 
in orthogonal directions. We will in the remaining sections take T""^ to be of the perfect fluid form and 
to be vanishing but we note that this worldvolume effective theory is valid for any type of material 
one would like to describe. 



2.2 Intrinsic fluid dynamics 

Our case study is that of a perfect fluid with energy density e{a°-) and pressure P{a^) with associated 
local entropy s^a"") and corresponding local temperature T{(t'^). The stress-energy tensor is given by 



nafe 



(e + P)n"n^ + P7 



Here the fields denote the fluid velocities normalized such that u°'Ua 
dynamic equilibrium the first law of thermodynamics must be satisfied 

de = Tds , 



(2.9) 

T. Assuming local thermo- 
(2.10) 



where the infinitesimal differentials are taken along worldvolume directions. In addition these quantities 
are supplemented by the Gibbs-Duhem relations 



e + P = Ts , dP = sdT . 



(2.11) 



The fluid dynamical equations follow directly from the intrinsic equation (2.7), which upon contraction 



along directions tangential and orthogonal to the fluid flows and using the above local thermodynamic 
relations can be written as the conservation of the entropy current J", 



DaJ! = 0, 



su 



and the Euler force equations 



r^Tsiiib + dblnT) = 



(2.12) 



(2.13) 



•;■ ab 



where _L = 7°^ + u'^u". Eqs. (|2.12|)-(|2.13l) are subject to the boundary condition (|2.6|) now in the form: 



,a„.b 



rj^ab^l. 



af^b\aw, 



P+l 



. 



(2.14) 
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Focusing on stationary fluid configurations, these two equations are solved by requiring the fluid velocities 
u"" to be proportional to a worldvolume Killing vector field k of the form |24| 

k = e + ^^^X^, (2.15) 

where ^ is a worldvolume timelike Killing vector, Xi ^-^^s spacelike ones and i labels spatial worldvolume 
directions. Without loss of generality we choose 

k" 

= — . (2.16) 

Jv 



The global fluid temperature T appears as an integration constant from Eq. (2.13) and it is related to 
the local temperature T via a local redshift factor, 

r = |k|r, (2.17) 

while the total entropy, assuming k^ to be hypersurface orthogonal with respect to the worldvolume 
metric, can be obtained by integrating the entropy current over the spatial part of the worldvolume Bp, 



dV^p^su^-Ua . (2.18) 

Here we introduced the spatial measure dV^p-^ on the worldvolume and defined an orthogonal vector Ua 
to a worldvolume spacelike hypersurface in the manner 

n» = 1^ , (2.19) 

where Rq is the norm of the timelike Killing vector field ^ on the worldvolume. 



2.3 Extrinsic elastic dynamics 



In this section we analyze the extrinsic dynamics of fluid configurations described by Eq. (2.8). We 
begin by defining the state of strain of the brane and the strain tensor and then proceed to describe the 
equations of mechanical equilibrium. 



The state of strain and the strain tensor 

The fluid configuration studied in the previous section lives on an infinitely thin surface described by 
the induced metric j^,^. We are interested in examinaning the thermodynamic properties of such a fiuid 
when the geometry is deformed along orthogonal directions. The metric 7^jy measures distances between 
neighboring points on the embedding surface, therefore, working under the assumption that the surface is 
thin and hence that variations in distances measured with _L^,^ can be ignorecj^ the metric 7^1^ describes 
the local state of strain of the brane. 

Let us define the state of strain of the brane prior to a deformation (unstrained state) by j^^. The 
length of infinitesimal spacetime distances along the surface is thus of the form 

ds'^ = j^^dxf'dx'' . (2.20) 

'^This is the usual assumption of classical elasticity theory when considering deformations of thin membranes and small 
strains and stresses |16| . 
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After a deformation, the state of strain is no longer described by jf^u but instead by the actual value of 
7^jy, hence the length of the infinitesimal element is changed to 

ds"^ = j^udxf'dx" . (2.21) 

Assuming the strains and stresses involved to be small, the difference between the length of the line 
elements of the strained and unstrained case along an arbitrary orthogonal vector Np is given by: 

As^ = ds^ - ds^ = £Nl,„.dx''dx'' . (2.22) 

Since we are only interested in changes along world volume directions we project this measure along those 
directions, 

As'lWp+i = lllP£NlXpdx^dx'' = -2NpKf,,Pdx^'dx'' , (2.23) 

where we have used a mathematical identity described in [lOj. Therefore, along any orthogonal direction 
Np, the strain varies proportionally to the extrinsic curvature tensor Kp,^'^ Since Kpi,^ satisfies the 
property Kp,yP = u'^pv!l,Kah'\ we define the Lagrangian strain tensor [25] for the brane as 

Uab = -\{lab-lab) , (2.24) 

which for infinitesimal deformations reads 

dUah = -\diab = NpKab" • (2.25) 
Equations of extrinsic dynamics 

The extrinsic dynamics of the fluid living on the dynamical surface are described by the extrinsic equation 



(2.8) which can be written in the form [10] 

Ts±Ppu'' + PRP = , (2.26) 

relating the acceleration of the fluid along orthogonal directions to the action of the mean extrinsic 
curvature Rf = j'^^Kab^- Here, ii^ is the fluid acceleration defined as = u^S/yU^. We now make 
the assumption of the existence of a background Killing vector whose pullback onto the worldvolume 
coincides with the worldvolume Killing vector field k". Under this assumption we can write 

Tsiip = -sdpT . (2.27) 



Introducing this into Eq. (2.26) and using the definition of RP yields. 



±PpsdPT = P^'KabP . (2.28) 



Contracting the last equation with an arbitrary orthogonal vector Np leads to 



sdT = -^a^^d^ab , (2.29) 



^Preliminary arguments in this direction were first given in |17) . 
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where the infinitesimal differentials denote a variation along orthogonal directions and where we have 
defined the pressure tensor 0-"'' as: 



(2.30) 



The r.h.s. of Eq. (2.29) can be written as a"'^dUab by recognizing the strain tensor defined in ( 2.24[ ). We 
now define the elastic solid density p as 



p = e + P = Ts 



(2.31) 



and hence interpret s as an average particle density satisfying the conservation law (2.12) and T as a 



mass function. In terms of infinitesimal variations along orthogonal directions to the worlvolume we find 



dp = Tds + sdT 

= Tds + a^^dUab = Tds - PdV , 



(2.32) 



where we have used the mathematical identity obtained in [10] for the relative change in the local volume 
element 



dV 



1 



(2.33) 



Eq. (2.32) is exactly the relation that an elastic solid density p should respect under hydrostatic com- 
and expresses the fact that the system described by this set of equations accounts for only 



pression 

changes in volume but not in shape. This is due to the assumption of an infinitely thin surface since 
effects due to bending or torsion would require a varying concentration of material across transverse 
directions. Along orthogonal directions we can thus write the equivalent relations 



dp 
dUab 



a 



ab 



dp 
dV 



-P . 



(2.34) 



Here, as well as in (2.32), we have assumed that all strain components dUat are linearly independent. 



This is not necessarily the case and we will deal with linear dependence towards the end of this section. 



We note that the first equality in (2.34) only constrains the components of the pressure tensor which 
are contracted with non-vanishing components of the strain tensor in (2.29). In order to make further 



contact with relativistic elasticity theory [251 [261 [27] note that by means of the definition (2.31 ) the above 
thermodynamic quantities and stresses can be derived from the mass function T, for example 



dT 
dV 



-P . 



Furthermore, we rewrite the stress-energy tensor (2.9) as 



(2.35) 



(2.36) 



motivating the interpretation of a solid at rest which has suffered hydrostatic compression along all 



worldvolume directions. The definition and interpretation of Eq. (2.31 ) together with the relations (2.32) 



and the decomposition (2.36) is one of the central results in this note as they express the elastic character 
of the relativistic material since the components of the pressure tensor a""^ involved in the extrinsic 
dynamics (2.29) can be obtained from a single potential p - a required condition in relativistic elasticity 
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theory [281 ESI EOj ETJ [30]. It is possible to consider a more general elastic potential which is also better 
suited when we consider charged branes. This will be the aim of the next section. For now, we focus on 



the case for which the strain components in (2.32) can be linearly dependent. In such situations, given 
the independent components of the strain tensor C/-^ we write (2.32) as 



dp = Tds + a~^^dU~^i , 
where we have introduced the effective pressure tensor along the component (a, 6), 

d^ab 



~ab I ab ab 



I db 
-L ab 



(2.37) 



(2.38) 



ab 



Here the operator -L"^a6 acts as a projector onto the linearly independent subspace of the components of 



strain. Using Eqs. ( 2.37 )-( 2.38) we can write 



dp 



ab 



a 



ab 



(2.39) 



We wish to rewrite an effective form for the stress-energy tensor (2.36) along the directions (a, 6). To 



this aim we note that from Eqs. ( 2.26 )-( 2.27) together with the identity u^u'^K^y'' = 1-'' ^ii [10] and the 
Gibbs-Duhem relation (2.11) one finds the expression 



pu'^u'' = 2 ' 

ab 



(97, 



(2.40) 



which is only valid along transverse directions. Acting with the projector -L"*a6 on both sides of (2.40) 



and using (2.38) we can write the effective stress-energy tensor along transverse directions as 



dP 

9^db 



~db 



which must satisfy the constraint: 



rpdb Q 



(2.41) 



(2.42) 



Imposing (2.42) leads directly to the extrinsic equations of motion (2.26). 



2.4 Elastic free energy, elasticity tensor and conserved charges 



Using the identity (2.11) it is possible to rewrite Eq. (2.29) only in terms of the pressure P as 

1 



dP 



-y'd^ab 



or alternatively, along orthogonal directions and for the independent components of d^y^^, 

dP 



9^db 



a 



ab 



(2.43) 



(2.44) 
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One can view this equation as a balance of forces between the pressure tensor and the internal stresses 
generated by a variation in volum^ Moreover, at the level of uncharged fluid branes, Eq. (2.44) is 
equivalent to Eq. (2.39) but this is not so in general. In fact, the pressure P, for reasons that will become 



apparent, provides a more general elastic potential which from now on we take to be the canonical one. 



From Eq. (2.43) it is possible to define the bulk modulus or modulus of hydrostatic compression K that 



measures the material response to variations in volume through the relation 

K ~ [dp)^ ~ ~P ■ 



(2.45) 



The definition (2.45) has a direct classical analog |16| . Eq. (2.43) can be integrated to an action [10 



I\X^'] 



(2.46) 



which resembles the usual action for fluids living on a fixed background, the difference being that 7 is 
the volume measure on the worldvolume instead of being on the space-time |31| . On the other hand, the 
worldvolume being a surface of co-dimension higher than zero, allows for elastic behavior which can be 
seen from the variation of the integrand along orthogonal directions 



-f'dP-d{^/^) P 
7 {-sdT-PdV) , 




(2.47) 



where we have used the Gibbs-Duhem relation ( 2.11[ ). Eq. (2.47) has direct analogy with the variational 
form of the Helmholtz free energy of an elastic solid, exhibiting the same local thermodynamic properties 
|16j . The action (2.46) demands to be interpreted as the solid free energy = — /[X'^] when the 



material has suffered hydrostatic compression. To motivate this interpretation even further we assume 
in the present moment the fluid to be barotropic, characterized by an equation of state e = wP. This 



together with (2.11) and (2.31) allows us to rewrite the solid free energy as 



F[X''] 



w + l 



(2.48) 



which is the usual action for relativistic elastic media |321 [261, l27] . The assumption of barotropy can be 



relaxed in order to obtain an expression of the form of (2.48) as we show in Sec. 2.6 when considering 
charged fluids. We note that the action (2.46) is very general extending also to charged fluid branes 
[33l [m [m [Ml [35] , while the action (2.48) is dependent on the equation of state. For this reason we have 



chosen the elastic potential P to be the canonical one instead of p. 



The elasticity tensor 

In elasticity theory, given the potential P from which the pressure tensor cr'^^ can be obtained, it is 
possible to define a deformation tensor K"'^'^'^ as the variation of the pressure tensor with respect to the 
state of strain [251 [26l [27] . For the fluid branes considered here this has the general form 



K 



abed 



da 



d'^p 



(2.49) 



^For the fluid branes arising from a gravitational dual analyzed in [TU] the pressure tensor takes the interpretation of 
gravitational tension acting as a compression force while the first term in Eq. (2.441 takes the interpretation of a centripetal 



force acting outwards the surface when the fluid is rotating. This is clear from the r.h.s. of (2.401 since it is proportional to 
a density p and two copies of the four- velocity 
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satisfying the properties K 



dbcd 



f^{^b){cd) ^ ^cddb_ ^^^^ ^i^g^^ ^j^g second equality in dHol) we have 



imposed the constraint (2.44). Given a certain stationary fluid brane in a certain prestrained state 



satisfying (|2.44[ ), one can apply a deformation taking the configuration to another strained state. The 
definition (2.49) implies that variations of the pressure tensor along orthogonal directions are related to 



the deformation tensor as 



da 



~db 



X^bcd^jj^ 



cd 



(2.50) 



In order to give a concise explicit expression for (2.49) we introduce an effective pressure P(7^^) along a 
direction (a, b) such that 



ab\^db 



Using then the definition (2.49) we obtain the general expression 



K' 



abed 




p^d{c^d)b 



(2.51) 



(2.52) 



This has the expected form of the deformation tensor of a material that is responding to stretching 
or compression signaling the fact that since the surface is taken to be infinitely thin, only variations in 



volume can be accounted for. The second term in (2.52 ) is the usual term when the pressure P is constant 
while the first term arises due to pressure variations. For Dirac branes the first term on the r.h.s. of 
(2.52) vanishes expressing isotropic compression at constant pressure. 

The relativistic elasticity tensor E°'bcd expressed in an analogous way with respect to the 

stress-energy tensor (2.41) \2b \ I26 t [27]. In orthogonal directions to the worldvolume, the elasticity tensor 
describes variations of T"*^ such that 



^fdb ^ E^bcd^jj__ ^ 



According to (2.53) the elasticity tensor takes the following generic form 



E 



dbcd 



dT 



'db 



dU,d. 



K 



dbcd 



(2.53) 



(2.54) 



where in the second equality we have used (2.52). The tensor E"'^'^'^ satisfies the usual properties of an 
elasticity tensor E"^''^ = = £;cJa6_ -y^g ^^^^ ^YlsX the Eqs. (|2.50l), (|2.53l) express hnear Hokean 



deformations of the pressure and stress-energy tensors. We will evaluate (2.54) explicitly for neutral black 
branes in Sec. ID 



Conserved charges 



The action (2.46) has a thermodynamic interpretation |10j . To make it precise, we note that from 



Eqs. (2.5 ) with vanishing external force one can construct a set of conserved surface currents T^'^k^, such 
that 



7^Vp(r^-k^) = 



(2.55) 
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where should be interpreted here as a generic space-time Kilhng vector field. The total energy M and 
angular momentum J* of the fluid brane can then be computed by integrating the surface currents over 
the spatial part of the worldvolume in the following way: 



M 



(2.56) 



Here we have assumed to be hypersurface orthogonal with respect to the space-time metric and also k^ 
to be hypersurface orthogonal with respect to the worldvolume metricQ To proceed further, we introduce 
the fluid Gibbs free energy density G as 



g = e-Ts = -P , 

which has the following thermodynamic properties along orthogonal directions 

dg = -sdT = PdV . 



(2.57) 



(2.58) 



From Eq. (2.58) we conclude that deformations of Q along orthogonal directions cause the material to 



stretch or compress. After integrating the density (2.57) over the worldvolume one finds the relation 



I ^Q = -\M-y^^,f-Ts\ . 



(2.59) 



to be satisfied for the fluid branes: 



Extremizing (|2.59|) while keeping the set of potentials T, fixed implies the first law of thermodynamics 

(2.60) 



dM = ^idf + TdS . 



This formula can interpreted as a prediction [TO], namely, that a stationary fluid configuration living on 
a particular dynamical surface must globally satisfy the first law of thermodynamics. 

2.5 Deformations of fluid branes and elastic waves 



The considerations of the previous sections are based on the action (2.46) being only a function of the 



state of strain given a fixed temperature T and angular velocities fij. Therefore, providing the pressure 
P as a function of the independent components of the state of strain ^ab-, i-e., P^'jab), deformations of 
the embedding geometry by explicit variation of the embedding map X^(cr") can be analyzed through 
variations of the induced metric jab- These take the form [33j 



Hence, searching for the extrema of the action (2.46) implies that 

T^'^Tab = , 

where the effective stress-energy tensor is obtained in the usual way 



lab 



61 



-7 H 



ab 



(2.61) 



(2.62) 



(2.63) 



*Other cases where or ka are not hypersurface orthogonal have been considered in |36l I35| . 
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For the present case, using (2.63), the stress-energy tensor has the general form 



-lab 



^ 

OJab 



(2.64) 



and is generahy different from a . Note that (2.64) is equal to (2.36) when the relation (2.40) is used 



The stress-tensor (2.64) obtained from the action is only valid along orthogonal directions and can be 



regarded as an off-shell form of (2.36) since a priori one can compute (2.64) without the knowledge of the 



extrinsic curvature of the embedding. The equations of motion obtained from Eq. (2.62) when projected 



along tangential and orthogonal directions to the worldvolume give rise to the intrinsic and extrinsic 



equations (1.1) [33] together with the boundary term 



l9w. 



(2.65) 



yielding the boundary condition (2.6) and hence vanishing independently of the initial and final configu- 



rations by construction. It is also possible to define an off-shell form of the elasticity tensor (2.54) using 

dP 



(2.64), this is given by 



abed 



2 P7"^'=7 



a{c d)b 



d-fabj yd^abdjcd 



(2.66) 



Projecting (2.66) using (2.38) for a particular embedding surface results in (2.54). 



The speed of elastic waves 

The dynamical properties of fluid branes as elastic materials can also be seen by applying a small pertur- 
bation to the brane geometry. To this aim we follow |10j . Assuming the material to be initially at rest 
= (1,0,...), we introduce a perturbation in the mapping functions 5X^ and initial pressure P such 
that 



5P , 5p= ^6P, du" = {0,v'), 5X^ = e 



(2.67) 



Using the form of the stress tensor (2.9) together with the equation of extrinsic dynamics (2.8) we find 

(2.68) 



[{p-P)d^^-Kd^) e = 

and hence conclude that elastic waves propagate at the speed 

K 



cl 



P 



P ' 



(2.69) 



where we used the definitions of the modulus of rigidity (2.45) and of the solid density (2.31). The above 



result is what is expected from a relativistic solid which has been subject to hydrostatic compression 

[371 ISg. 
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2.6 Charged fluid branes 



Fluid configurations carrying a g-charge are not only characterized by the stress tensor (2.3) but also by 



a totally anti-symmetric current tensor jMi-'-Mg+i which is confined to the worldvolume surface 

J^i-^«+i(x°) = / r/P+i^^/=YJ^i-^^+i(^°) " ^ >> - (2.70) 



and must satisfy current conservation 



l = , (2.71) 



in the absence of external couplings. The set of effective worldvolume equations that result from (2.71) 
imply that the current is purely tangential 

and 

jj^ja,...a,+, ^ Q ^ jai...a,+i^^j^^^^^ ^ g _ (2.73) 

We will now apply these equations to the fluid configurations with q = p charge studied in |13) and with 
q = charge studied in [TH [13] . 

q = p worldvolume charge 

The fluids carrying q = p worldvolume charge studied in pjij are characterized by a worldvolume current 
density J'^i- '^p+i of the form 

J = QpVp+i , (2.74) 
where V^+i is the {p + l)-volume form of the embedding surface and Qp is the charge density. The 



worldvolume conservation equation (2.73) then implies 

daQp+i = . (2.75) 

Therefore for q = p charges the charge density Qp is not allowed to vary along worldvolume directions, 
hence the total charge of the configuration Qp equals the charge density Qp. Given this, the fiuid does not 
carry any additional extra degrees of freedom associated with the charge and as such the Gibbs-Duhem 



relations presented in (2.11) still hold for the case at hand. In this way, as long as all extrinsic variations 
are performed while the charge Qp is kept constant, the results of the previous sections hold. However, it 
is possible to introduce a chemical potential $p conjugate to Qp [13], giving rise to a well defined global 
quantity 

= - I dV(p)$pn'^n, . (2.76) 

Using $p we can define the Gibbs free energy of the fluid as 

g = e-Ts- %Qp = -P- %Qp , (2.77) 
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where in the second equahty we have made use of the relations (2.11). Therefore, for hydrodynamical 
fluctuations where the charge Qp is kept constant due to (2.75) we find 



dg 



-sdT - Qjyd^ 



p ' 



(2.78) 



and hence for orthogonal variations that keep the potential (2.76) constant, ie., 

(i$p = - f^pdV , 



(2.79) 



the action can be recast as (2.59) noting that along those directions dQ = —QdV. Prom here one can 
define an electroelastic modulus of rigidity describing the deformation of the charge potential as 

av \ _ _^ 
)t~ ^p' 



1 



(2.80) 



As far as the authors are aware, Eq. (2.80) does not have a classical analog. We suspect that this modulus 
of electroelastic rigidity is associated with fluctuations of the charge density in transverse directions in 
the same way as the isothermal permittivity is associated with worldvolume fluctuations [13] . We further 



note that it is possible to define an electroelasticity tensor associated with $p as in (2.54) but its meaning 
is unclear. A better understanding of this is presently lacking. 

Finally, extremizing (2.59) at fixed T, <I>|^^ implies the first law of thermodynamics |13j 



dM 



VLidf + TdS + ^^^UQp 



(2.811 



to be satisfied. 



g = worldvolume charge 

Stationary fluids carrying a g = brane charge were analyzed in |14[ [T3] and have quite different 
thermodynamic properties than the q = p case. These are instead characterized by the worldvolume 
particle currenlj^ 



and hence must satisfy (2.73): 



(2.82) 



(2.83) 



The crucial difference with the q = p case is that now the charge density Q is allowed to vary along the 
worldvolume and hence adds extra degrees of freedom to the system. Local thermodynamic equilibrium 
implies 



de = Tds + 



while the thermodynamic Gibbs-Duhem relations (2.11) are now changed to 

e + P = Ts + <^Q , dP = sdT + <^>dQ . 

^From hereon we omit the index in Qp and "I>p. 



(2.84) 



(2.85) 
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The intrinsic equations of motion (2.7) again lead to conservation of the entropy current as in (2.12) 
while the Euler equations (2.13) are modified to [14j 

P^^Tsiiib + dblnT) - Q$ - P^^db In = , (2.86) 

where we have introduced the mean curvature of the worldlines embedded in Wp+i, 

i^'^ = uVDfe (ucu") . (2.87) 

Assuming the solution to be stationary, we take the fluid velocities to be aligned with a worldvolume 
Killing vector field as in (2.16). This ensures that the first term in Eq. (2.86) vanishes leading to the 
same relation between local and global fluid temperatures (2.17). The second term in Eq. (2.86) can be 
dealt with by imposing the mean curvature due to the dissolved 0-charge on the worldvolume to balance 
the gradient of the chemical potential [Hj , yielding 



K" = P''''db\n^ 



On the other hand, the extrinsic equation (2.27) now becomes 

Ts±P^u'' + Q±%d^<^> + PRP = , 



which after using (2.16) and contracting with an arbitrary orthogonal vector Np leads to 

PdV + sdT + Qd^ = . 



(2.88) 
(2.89) 
(2.90) 



This can be integrated to the action (2.46). In order to make further contact with electroelasticity we 
define the solid density 



which under deformations along the extrinsic directions satisfies 

dp = Tds - PdV + ^dQ . 



(2.91) 



(2.92) 



This is the thermodynamic relation that a solid charged under a particle current should respect. From 
here, as before, we can obtain useful relations along orthogonal directions: 



dp 



(2.93) 



Unlike the neutral case, not everything can be derived from a mass function. In fact, if we use the 



definition (2.91) into (2.93) we find the identities 

'or 

dV 



+ Q 



df 
[dV 



s 

'Q ' 
-P . 



(2.94) 



The Gibbs free energy introduced in (2.77) for charged fluids is now, due to (2.85), equal to the pressure 

(2.95) 



as in (2.57). Hence, if one wishes to specify an equation of state of the form 

e = wP + , 
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one finds the action (2.48): 



1 



w + l 



(2.96) 



Wp+i 



When written in terms of the sohd density p, the extra term on the r.h.s. of the action (2.96) acts as an 
external force. Variation of (2.96) keeping T, J7i,<I>H constant leads to the first law of thermodynamics 



(2.81), but now with global charge ^14j : 



Q 



(2.97) 



3 Application to stationary black j9-branes 

In this section we apply the considerations of the previous sections, as an example, to the case of neutral 
black p-branes wrapped on a generic submanifold (blackfolds) studied in [HI [TUl ESI EHl HO] . Such branes 



are characterized by a stress-energy tensor of the type presented in (2.3) and hence encompassed within 



the framework put forth above. We begin by describing the effective blackfold fluid characterizing these 
branes and then write down the elasticity tensor as well as an action of the thermodynamic type for 
blackfold objects. We finish by re-deriving the elastic equilibrium condition for higher dimensional black 
rings and black odd-spheres 



3.1 The effective blackfold fluid 



The blackfold effective fluid characterizing neutral black branes is of the perfect fluid type (2.9) and 
characterized by the equation of state [TO] 



-{n + l)P 



(3.1) 



The local thermodynamic quantities associated with the fluid living on the brane (e, P, T, s) can all be 
described in terms of the brane thickness tq in the following way: 



n + l 



-P 



(n+l)^0 

167rG 



T 



n 



47rro 



o n+l 
"(n+l)^0 

4G 



(3.2) 



Here n is the number of transverse directions given by n = D — p — 3. By inspection of the above 



expressions it is easily observed that the Gibbs-Duhem relations (2.11) are satisfled. The thickness tq 



as well as the fluid velocities u"' are allowed to be functions of a" and hence to vary along worldvolume 



directions. Prom (3.2) one concludes that the requirement of local thermodynamic equilibrium (2.1) 



implies the hierarchy of scales 



ro(a'^) « i?(a«) 



(3.3) 



or in other words, that the black branes being wrapped on the submanifold Wp+i are thin compared to 
the curvature radius of the submanifold. 



16 



3.2 Elasticity tensor of blackfolds and thermodynamic action 



Here we write down the elasticity tensor for neutral blackfolds using Eq. (2.66) and then an action of the 



thermodynamic type. We begin by noting that from Eq. (2.17) and using the quantities (3.2) we obtain 
a relation between |k| and rg 



n 



A|k| , 



(3.4) 



where we have defined A = n/ATiT. Using this, the pressure P can be expressed in terms of the induced 
metric as 



167rG 

Hence, the derivative of the pressure with respect to the induced metric is simply given by 



97, 



ah 



Using the expression for the effective stress-energy tensor derived in (2.64) we arrive at 

rpab 



P 



(3.5) 
(3.6) 
(3.7) 



agreeing with (2.9) when the equation of state (3.1) is introduced. In order to continue further, it is 

(3.8) 



necessary to evaluate the second derivative of P with respect to the state of strain. This yields: 

a^P \ n(n-2 



cd 



-Pu^vl'u^'^ 



Therefore, using Eq. (2.66) we obtain the off-shell elasticity tensor for neutral blackfolds: 

P {2^<^-f'^)^ + nu^u^-f"^ - n{n - 2)u^u''u^u'^^ . 



E' 



abed 



(3.9) 



This is not manifestly invariant under (a, 6) — >• (c, d) as discussed around Eq. (2.54) but it becomes so 



after the equations of motion (2.44) are imposed. We note that (3.9) has the same structure as the Young 
modulus tensor measured for black branes in |17l \V2\ lisj^ Moreover, it satisfies the properties 



jpabcd.^ 
-C/ 7cd 



P ( 27"^ + n(n + p - l)'u"u' 



Tpabcd^, „ 
-C/ gabled 



{2{p+l) -n{n+p-l))P . (3.10) 



Thermodynamic action 



Here we rewrite the action (2.46), adapted to blackfold objects, in terms of the total entropy (2.18). 



Using the Gibbs-Duhem relations (2.11) and also (3.1) we obtain the relation 

P = -Ws , 



n 



which when introduced into the action (|2.46|) leads to 



— / dV(„)SU Ha = 



(3.11) 



(3.12) 



Therefore we see that extremizing the blackfold action at constant temperature and angular velocities is 
to extremize the black hole entropy. 



There is a difference between the structure of (3.91 and the Young modulus measured in |17l 1121 [TS] . namely that here 
there is no ambiguity in the choice of worldvolume surface since we take the surface to be infinitely thin. 
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3.3 Equilibrium condition for black holes 



In this section we apply the results of the previous sections to re-derive the equilibrium condition of 
thin black rings and black odd 3-spheres [38] in asymptotically flat space parametrized by coordinates 
(t, r, 6, (j), ^, x^). We write the induced metric on the worldvolume as 



2 + cos^^d^^ + sin^^dV^) , 



(3.13) 



where R is the radius of the sphere (or the ring). The map onto the ambient space-time is simply 

t = T, r = R, x' = , (3.14) 
while the remaining coordinates coincide. The non- vanishing components of the extrinsic curvature are 

Kee'' = -R K^/ = -Rcos'^9, K^/ = -Rsin^O , (3.15) 
while the mean extrinsic curvature vector reads 



R 



p=l,3 



(3.16) 



Black rings 



Black rings are described by setting p = 1 and 9 = in the formulae above. This leads to the induced 
metric and extrinsic curvature 



-fabda'^da'' = -dr'^ + R 



2 Jj:,2 



-R , 



1 

R 



The worldvolume Killing vector field takes the form 

k = dr + nd, 

which we use to write the pressure P as 



p 



levrG 



(3.17) 



(3.18) 



(3.19) 



Since the only independent and non-zero component of the strain tensor is 
obtain the equilibrium 



we can use Eq.(2.44) to 



dP 



0, 



which in turn implies, using the result (3.6), 



n+l 



(3.20) 



(3.21) 



This equilibrium condition has been derived previously in [38] also using the action (3.12), here we have 



merely performed a different derivation. From Eq. (3.20) we see that there is only one non- vanishing 



component of strain, hence the only on-shell component of the elasticity tensor is given by (3.9): 



(^) 



(3.22) 
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Black odd-spheres 

Black odd 3-spheres are described by setting p = 3. We assume them to be rotating with equal angular 
velocity ^1 along both directions {(p,ip). This allows us to write the worldvolume Killing vector field in 
the form 



and hence the pressure as 



IQttG 



A" 



7tt - iee^^\ 2 



There are now three non- vanishing components of the strain, d'yee, 
expressed in terms of d'yeo in the following way: 

c^700 = cos^Od-yee , d'y^^ = sin^ 9d'jg 



(3.23) 
(3.24) 

and dj^^. The last two can be 



(3.25) 



Using now Eq. (2.44) for the only independent component d'ygg we obtain the equilibrium equation 

dP 



d-fg 



where the effective stress along djgg is given by 



a 



P + cos^ 



a 



+ sin" 







Solving Eq.(3.26) results in the equilibrium condition 



p 

n + p 



(3.26) 



(3.27) 



(3.28) 



which has been derived previously in [38]. Even though we only considered the cases p = 1,3, the result 
(3.28) is valid for all p [38]. Eq. (2.44) could have been solved using any of the other two non- vanishing 
components of strain U^^,U^jp, in fact, it is easy to see that the effective pressure tensor has in this 
case the generic form a""^ = P^°-^ where P = pP. Using expression (2.54) for the on-shell value of the 
elasticity tensor we find the non-vanishing components: 



E' 



n + p 



n 



n + p 



n 



2 ( ] p^^i'^^i' 



n 



2^P7^S<^ 
n 



E' 



n 



2— p^<t>(t>^-4'i' 



(3.29) 



n 



4 Conclusions 

Relativistic elasticity theory has been considered in the literature by a number of authors [281 132| W\\ [29l 
W2\ [26] [27] , one of the most influential works being that conducted by Carter and Quintana [25] . The 
main focus of these works has been on describing elastic solids which are space-filling, contrary to being 
confined to a dynamical surface. However, the main difference compared to the work presented here is 
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not only related to the existence of an embedding surface but also to the usual assumption regarding the 
form of the pressure (2.30) and elasticity (2.54) tensors: 



a''^ = , E'^^^'^Ud = . (4.1) 
This orthogonality condition with respect to the fluid flows seems to have classical roots, namely, that an 



elastic material can be deformed in space but not in time. Not imposing (4.1) has been the main point 
of departure in this note from the usual formulations of relativistic elasticity and has drawn inspiration 
from the recent results of bending black branes \17\ [T2| [T8] . 



Prom the point of view of material science, materials that do not satisfy (4.1) might belong to a 
rather unusual class, but from the point of gravitational physics, in which black holes in certain regimes 
may take the effective description as a fluid brane [6]- [12], they are relevant for the analysis of their 
dynamics. Their study can lead to the formulation of a corner of relativistic elasticity theory that has 
not been previously explored (e.g. by applying extrinsic perturbations to black brane geometries) in the 
same way that the fluid/gravity led to a more general formulation of dissipative corrections to fluid and 
superfluid hydrodynamics EOj EU [191 Ell ES]. Here, as an example, we have written down the most 



general form of the elasticity tensor of a fluid brane that responds only to changes in volume (2.54) and 



computed it for the particular case of neutral blackfolds (3.9). The assumptions taken were that the 
intrinsic properties of the material living on the thin surface were those of a stationary fluid in local 
thermodynamic equilibrium whose backreaction onto the ambient space-time is negligible. It is unclear, 
at the present moment, wether taking another kind of intrinsic dynamics would lead to an elastic behavior 



since it was crucial for our understanding that the extrinsic equations of motion (2.26) could be integrated 
into an action that only depends on the distances between neighboring points on the brane after fixing 
the potentials T, Qi, <I>h- 

It is not surprising that a stationary fluid will behave elastically along orthogonal directions to the 
surface on which it lives. Given a perfect fluid with density e and pressure P living on a two-dimensional 
sphere embedded in four-dimensional space-time and applying a deformation to the surface such that the 
radius of the sphere increases by a small amount, then the density and pressure will drop (or increase) due 
to the increase in volume. This is what we observe here. It will be interesting to apply these ideas and 
develop the theory when finite thickness corrections to the surface are considered \n\ [T8] . Considering 
such effects allows for probing other types of deformations, such as bending, which require a varying 
concentration of matter along transverse directions. 
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